Critical properties of the Fermi-Bose Kondo and pseudogap Kondo models: 

Renormalized perturbation theory 
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Magnetic impurities coupled to both lermionic and bosonic baths or to a fermionic bath with 
pseudogap density ol states, described by the Fermi-Bose Kondo and pseudogap Kondo models, dis- 
play non-trivial intermediate coupling fixed points associated with critical local-moment fluctuations 
and local non-Fermi liquid behavior. Based on renormalization group together with a renormalized 
perturbation expansion around the free-impurity limit, we calculate various impurity properties in 
the vicinity of those intermediate-coupling fixed points. In particular, we compute the conduction 
electron T matrix, the impurity susceptibility, and the residual impurity entropy, and relate our 
findings to certain scenarios of local quantum criticality in strongly correlated lattice models. 



I. INTRODUCTION 

Zero-temperature phase transitions in quantum impu- 
rity models have attracted significant interest in recent 
years in a number of different contexts. In correlated 
electron materials impurities have proven useful for in- 
vestigating the complex bulk behavior, by probing the re- 
sponse to point defects e.g. in a spatially resolved manner 
by NMR or scanning-probe microscopy techniques. In 
quantum dot systems, impurity models naturally arise in 
describing the degrees of freedom on one or multiple dots 
together with their coupling to leads. Finally, the physics 
of strongly correlated bulk systems can be mapped to im- 
purity models in the framework of the dynamical mean- 
field theory (DMFT)A, where a local approximation to 
the lattice electron self-energy leads to an impurity model 
supplemented by a self-consistency condition. 

Interesting transitions, where the zero-temperature be- 
havior of the impurity degrees of freedom changes qual- 
itatively as a function of some control parameter like 
coupling constants or gate voltages, occur in modifica- 
tions of the familiar Kondo model, which describes an 
impurity spin interacting with conduction electrons of 
a metallic environment^. Whereas the standard Kondo 
model does not show a phase transition at a finite value 
of the coupling, transitions naturally occur when other 
effects, which compete with fermionic Kondo screening, 
are added. At the resulting transition point, the Kondo 
energy scale is suppressed to zero, and the low-energy 
physics is dominated by critical local-moment fluctua- 
tions leading to singular thermodynamic properties and 
local non-Fermi liquid behavior. Notably, the physics at 
these intermediate-coupling fixed points is different from 
the well-studied multi-channel Kondo behavior. 

In this paper, we will focus on two routes to impu- 
rity phase transitions. The first is represented by the 
so-called Fermi-Bose Kondo model'^-'*-^, describing an im- 
purity spin interacting with both fermionic quasiparti- 
cles and collective spin fluctuations. The Fermi-Bose 
Kondo model has recently been analyzed for the case of 
a metalHc fermion density of states (DOS), and a gap- 
less bosonic spectrum representing magnetic order pa- 



rameter fluctuations at a bulk quantum critical point in 
d dimensions^i^. For d < 3 the model shows an impu- 
rity quantum phase transition between a Kondo-screened 
phase and a bosonic fluctuating phase with universal lo- 
cal spin correlations. Recently, the Fermi-Bose Kondo 
model has been proposed^ to describe a "local quan- 
tum phase transition" in alloys like CeCug-ajAuj^. This 
modelling is based on an extended dynamical mean-field 
theorji^, where the impurity site is coupled to both a 
fermionic and a bosonic self-consistent bath. 

A second paradigmatic model is the pseudogap Kondo 
modelSi^iifliiLiSii^ where the fermionic bath displays a 
pseudogap DOS, pd'^) oc ■ Such a bath arises in 
particular in ordered fermionic systems where the order 
parameter has nodes at the Fermi surfaces, like unconven- 
tional superconductors. In the pseudogap Kondo model, 
screening is suppressed for small Kondo couplings, and a 
transition between a Kondo-screened and a free-moment 
phase results. The pseudogap Kondo model has been 
extensively studied by a number of techniques, in partic- 
ular using the numerical renormalization group method 
(NRG) . Recently, it has become clea ri . 4S that r = and 
r = 1 play the roles of the lower and upper-critical "di- 
mensions" , respectively, allowing for different perturba- 
tive expansions to access the physics of the critical point. 
On the application side, the pseudogap Kondo effect is 
relevant for impurities in d-wave high-temperature super- 
conductors, where indeed non-trivial Kondo- like behavior 
has been observed connected with the magnetic moments 
induced by Zn and Li impurities in YBa2Cu307_,fi^4i^ii&. 
Furthermore, it has been suggested-^'' that quantum dots 
coupled to interacting one-dimensional electron liquids 
can display pseudogap Kondo physics. 

It is apparent that an analytical description of the 
properties near the intermediate-coupling fixed points is 
highly desirable to obtain a better understanding of the 
physics of critically fiuctuating moments. In this paper 
we shall calculate a number of properties of the Fermi- 
Bose Kondo and the pseudogap Kondo models for param- 
eters in the vicinity of the intermediate-coupling fixed 
points. These calculations are based on a perturbative 
renormalization group analysis of the models, together 
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with renormalized perturbation theory - this allows us 
to obtain results for observables in the quantum-critical 
regime where bare perturbation is known to diverge. As 
in earlier work^i^iSiiSiiSiSS, the expansion is done around 
the free-moment fixed point (which can be understood 
as expansion around the lower-critical dimension of the 
phase transitions). The small parameters controlling the 
perturbation theory are the exponents r and e = 3 — c? 
occuring in the fermionic and bosonic bath density of 
states, respectively. Critical properties will be computed 
in a double expansion in the non-linear couplings of the 
theory and the scaling dimensions r and e. In partic- 
ular, we will show that the impurity entropy is finite at 
the Fermi-Bose critical point, which has consequences for 
scenarios of local criticality in the framework of extended 
DMFT. 

The paper is organized as follows. Sec. ^1 will intro- 
duce the impurity models of interest and describe the 
properties of the fermionic and bosonic baths. In Sec. lIIII 
we will briefly outline the renormalization group treat- 
ment which is based on a weak-coupling expansion of the 
impurity-bath couplings. In particular, all fixed points 
of the general pseudogap Fermi-Bose Kondo model will 
be enumerated and discussed. The calculation of observ- 
able impurity properties, namely susceptibility, entropy, 
and T matrix, is presented Sees. IIVIIVII When possible, 
the results are compared with numerical data. Finally, 
Sec. I VIII discusses implications of our study for scenar- 
ios of local criticality in heavy- fermion metals. Technical 
details are relegated to the appendices; in addition we 
present a large- TV calculation of the impurity entropy for 
the Bose Kondo model in App.IdI which is not restricted 
small e and thus complements the e expansion of Sec.lVl 

Parts of the RG analysis and susceptibility calculations 
have been presented in Refs. 1511 11912012 It here we general- 
ize these earlier works and derive additional results. The 
main new findings concern the impurity entropy, which, 
to our knowledge, has not been explicitely evaluated be- 
fore and is an interesting property in the DMFT context; 
in addition, we present magnetic response results for the 
pseudogap Kondo model which can be nicely compared 
to numerical data. 



II. MODELS 

It is straightforward to combine the effects of a pseu- 
dogap bath and a competing bosonic interaction onto 
fermionic Kondo physics by considering the pseudogap 
Fermi-Bose Kondo model, which supersedes the variants 
discussed above. The Hamiltonian of the quantum im- 
purity problem can be written as Ti. — Tic + Ti,^ + Himp, 
with 

Himp = 7o<S'Q0a(a; = O) + joS'c<(4(T^'J^^c^)(?; = 0) . (1) 

Here, S is the operator of the impurity spin, with 
[Sa,Sj3\ = iea/s-ySj, and we will restrict our attention 



to 5 = I in this paper. The fields (f>a{x) and c^{y) de- 
scribe bulk spin fluctuations and bulk fermionic excita- 
tions, respectively, with Hamiltonians Ti,^ and He to be 
discussed in the subsections below. In particular, the 
fermionic bath will follow a power law density of states, 

Pc(w) (X i^r , (2) 

which includes the cases of a metal (r = 0) and a d- 
wave superconductor {r — 1). v, fi are spin indices with 
ly =1, 1, cr(") is the vector of Pauli spin matrices. 

The impurity part of the Hamiltonian Tiimp, contains 
couplings to the bosonic and fermionic baths, character- 
ized by coupling constants 70 and jo. In the absence of 
the bosonic bath, 70 = 0, and for r ~ the model 
reduces to the well-known Kondo model-, which exhibits 
screening of the magnetic moment below a characteristic 
(Kondo) temperature, Tr-. The pseudogap generaliza- 
tion, r > 0, shows a zero-temperature phase transition 
between a free local moment and a screened moment at 
a critical coupling jo — joc- On the other hand, in the 
absence of the fermionic bath, jo — 0, the Bose-Kondo 
model describes an impurity moment in an insulating 
quantum antiferromagnelsiSiiSiSi; in addition it arises in 
mean-field theories for spin liquids and spin glasses^^. 
Our analysis will be restricted to a model with full rota- 
tion symmetry in spin space; the effect of spin anisotropy 
for r = has been discussed in Ref. |3 

A. Fermionic bath 

The fermionic bath consists of non-interacting spin-i 
fermions. (More precisely, the fermionic self-interaction 
is assumed to be irrelevant in the RG sense.) A conve- 
nient continuum representation uses chiral fermions with 
a linear dispersion in (1 + r) dimensions: 

nc = J dfc|A:r(«j.fc)cLcfc.. (3) 

The field c,y{y = 0) in Eq. ^ is simply the Fourier trans- 
form of Cfei,, Ci/(y = 0) = / <ik\k\^Cki^- The velocity vp 
will be set to unity in the following. 

We assume particle-hole symmetry of the bath; in the 
applicability range of our expansion particle-hole sym- 
metry breaking is marginal for r = and irrelevant for 
r > and thus does not change our conclusions. As will 
be mentioned in Sec. IIIII particle-hole symmetry break- 
ing becomes important for larger r and/or in the strong- 
coupling regime. 

In c?-wave superconductors, r = 1, and the c fermions 
represent the superconducting Bogoliubov quasiparticles 
in the vicinity of the nodes of the d-wave gap. 

B. Bosonic bath 

The bosonic bath, representing antiferromagnetic spin- 
1 collective fluctuations of the host material, will be de- 
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scribed by the 0(3)-synimetric quantum cj)'^ theory in d 
dimensions, 



neb = / d'^x 



2 +4! 



, (4) 



where the field (j)a (x) represents the local orientation and 
magnitude of the (non-conserved) antiferromagnetic or- 
der parameter, and its canonically conjugate momentum 
is TTa{x), hence [0ct(x), 7r^(a;')] = iSai36'^{x — x'). The 
parameter s tunes a quantum transition between a para- 
magnetic phase for s > Sc and an ordered phase breaking 
the 0(3) symmetry at s < Sc- 5o is the self-interaction 
of the </) quanta, and c is a velocity, and we will use units 
such that c = 1 in the following. The phase transition of 
the 0* theory can be studied by working in (3 — e) space 
dimensions, using a perturbative expansion in e. 

The quantum paramagnet for s > Sc is characterized 
by its spin gap Ag at T = 0, and this vanishes as s 
approaches Sc from above as 



As OC (s - Sc) 



(5) 



where i/0 is the correlation length exponent of the mag- 
netic bulk transition described by Ti.^; the dynamical ex- 
ponent is z = 1. The local two-point (j)a propagator in 
imaginary time r is given by 



D4t) = (0a(O,r)</)„(O,O))o 



T 



E 



d'^fc 



fc2 



(6) 



Here the momentum k is measured relative to the order- 
ing wavevector Q of the magnet. The subscript in the 
correlator indicates that it is evaluated to zeroth order in 
go- However, we have included Hartree-Fock renormal- 
izations in the determination of the "mass" m; these have 
been computed in earlier work in the e expansion^, and 
we quote some limiting cases to lowest order in e = 3 — d: 



A, 



v/10e/337rT 



s > Sc, T < As 
T> Is-ScT* 



(7) 



It is useful to define a spectral density, ptpioj), of the 
Bose field 0. At T = it will display a gap of size Ag; 
for As = it follows p,/,(a;) oc sgn(a;)|a;|-^~'^. 

For d < 3 the bosonic self-interaction g^ is relevant 
in the RG sense and cannot be neglected. However, in 
the RG analysis of the impurity model CQ), modifica- 
tions from go arise only at two-loop order. Treatments of 
the Fermi-Bose Kondo model in the context of extended 
DMFTi^ worked with a go = theory - there the irrel- 
evance of the bosonic self-interactions was argued to be 
justified due to the Landau damping of the spin fiuctu- 
ations. In what follows, we will treat both g^ — Q and 
go 7^ simultaneously, and display all results for both 
cases. In applications, the go 7^ case is relevant in the 
absence of Landau damping of the spin fluctuations, i.e., 
for metallic systems with a small Fermi surface obey- 
ing IQI > 2kp and for pseudogap systems like d-wave 
superconductors . 



III. RENORMALIZATION GROUP 

The Fermi-Bose Kondo model Q can be conveniently 
analyzed using standard renormalization group tech- 
niques. We shall assume that the spectrum of the bosonic 
bath is gapless. As = at T = 0, which corresponds 
to a magnetic bulk quantum critical point. Then, both 
fermionic and bosonic baths obey a power law DOS. The 
RG can be formulated as e xpansion in the bath exponents 
e and r, similar to Refs. I3l l5ll8ll9l 

It proves convenient to switch from the spin operators 
for the impurity in to an Abrikosov pseudo-fermion 
representation of the impurity spin ^, S — f},(Ti,fj,ffj,. 

The required Hilbert space constraint flfv ^ Q = ^ 
will be implemented using a chemical potential A ^ 00, 
such that observables {O) have to be calculated according 
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(O) 



lim 



(8) 



where (. . .)a denotes the thermal expectation value calcu- 
lated using pseudo-fermions in the presence of the chem- 
ical potential A. Clearly, in the limit A — > 00 the term 
{Q)\ represents the partition function of the physical sec- 
tor of the Hilbert space times exp(— A/3). As detailed in 
App.El both numerator and denominator of Eq. JHJ have 
to be expanded in the non-linear couplings to the required 
order when calculating observables; however, the denom- 
inator does typically not develop logarithmic singularities 
at the marginal dimension. 

As an alternative to the pseudo-fermion technique, the 
approach of Ref. can be utilized, which does not em- 
ploy pseudo-fermions, and is completely equivalent for 
the present single-impurity problem. 



A. RG transformation 

Some crucial properties of the partition function Z = 
Tr e"''^ follow from tree-level power-counting. The be- 
havior of the fields under a rescaling transformation de- 
termines the scaling dimensions: 



dim[0„] = (d-l)/2; dim[ci,] = (1 + r)/2 ; 
dim[/,] = 0. 



(9) 



The last line refers to the pseudo-fermion representation 
of the impurity spin. The above can be used to determine 
the dimensions of the couplings: 



dim [go] 
dim[7o] 



3- 
(3 



-d)/2- 



dim [jo] 



(10) 



Thus, the bosonic bulk interaction go is relevant for 
d < 3, i.e., the bulk critical point will be characterized by 
a finite fixed point value of g. The bosonic impurity cou- 
pling 7o is a relevant perturbation about the decoupled 
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FIG. 1: Schematic RG flow diagram for the pseudogap Fermi- 
Bose Kondo model for smaU e and r, and As = 0, deduced 
from the beta functions 1121 1. llAllll . The horizontal axis de- 
notes the coupling to the bosonic bath, 7, the vertical axis 
denotes the fermionic Kondo coupling, j. The solid dots are 
(meta)stable fixed points, the open dots denote critical fixed 
points, the thick line is the separatrix, for details see text. 



B. Fixed points 

In the following we list the RG fixed points, obtained 
from HA11|I . and quote some of their observable proper- 
ties like the Curie part of impurity susceptibility Tximp 
and the residual impurity entropy 5i,np; details of the 
calculation appear in Sees. llVl andlVl 



1. Local-moment fixed point (LM) 

The decoupled impurity fixed point is simply charac- 
terized by 7* = j* = 0. The susceptibility shows the full 
Curie moment of a spin i, Tximp = 1/4, and the impu- 
rity entropy is S'imp = ln2. 



impurity fixed point in d < 3 - this will yield a stable 
fixed point with finite T^iiSiiS. In contrast, the fermionic 
coupling jo is irrelevant, which eventually results in an 
unstable fixed point describing the phase transition in the 
pseudogap Kondo model. The combined effect of 70 and 
jo will be captured in the following RG treatment. 

To one-loop order, the familiar momentum shell 
method with ultra-violet (UV) cutoif A can be em- 
ployed; the field-theoretic RG calculation to two-loop or- 
der will be sketched in App. ^ Within the momentum 
shell scheme, dimensionless couplings are introduced as 
.90 = .9A75<i+i, Jo = JA^^ 70 = iK^'^S^'^^ with the 
phase space factors Sd, Sd defined in Eq. (|A6I) . Without 
loss of generality the same value of the UV cutoff can be 
used for the fermionic and bosonic subsystems. The RG 
for the bulk bosonic theory gives 



-eg- 
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(11) 



2. Fermionic strong-coupling fixed point (F-SC) 

For large fermionic Kondo coupling j, the bosonic part 
7 is irrelevant, and the model reduces to the fermionic 
pseudogap Kondo model which shows strong-coupling 
phases with screening (j* = 00, 7* = 0). This regime 
is not accessible within the present perturbation expan- 
sion, and we just quote the known results for the Kondo 
model. For the metallic case, r = 0, there is a line 
of strong-coupling fixed points parametrized by particle- 
hole asymmetry, along this line T^imp = and S'imp = 0. 
In contrast, for r > two different strong-coupling fixed 
points exis'feii. The particle-hole symmetric one is unsta- 
ble w.r.t. particle-hole symmetry breaking and only ac- 
cessible for r < 1/2, it does not display complete screen- 
ing: the susceptibility follows Tximp = ?'/8, and the im- 
purity entropy is S'imp — 2rln2. In contrast, the asym- 
metric strong-coupling fixed point is stable and exists for 



all r > 0, with Txm 



and Sii 



0. 



For the impurity part, the one- loop RG result for the 
beta functions is 



3. Bosonic fluctuating fixed point (B-FL) 



■2 , ■ 2 



(12) 



the two-loop results can be found in App.^in Eq. IjAlip . 
Already at one-loop level, a term mixing j and 7 appears, 
which determines the competition between fermionic and 
bosonic Kondo coupling. The resulting RG flow in the 
j-7 plane is shown in Fig. ^ the fixed points will be 
discussed in the following subsection. 

The equations H12|l . I|A11|I reduce to the cases of the 
pseudogap Kondo model for 7 = (Ref. the Bose- 
Kondo model for j — Q (Ref. 19), and the metallic Fermi- 
Bose Kondo model with non-interacting bosons for ?' = 
and 5 = (Refs.m. 



In the absence of a fermionic bath, the Bose Kondo 
problem shows interesting behavior: it flows to a stable 
intermediate-coupling fixed point, characterized by 



V 121 



132 



0(e3 



) -f 0{e^) for 50 ^ 
for 50 = 



= 0. 



(13) 



At this fixed point the fermionic Kondo coupling j is ir- 
relevant. The intermediate-coupling nature leads to uni- 
versal local-moment fluctuationsiSiiS, with T^imp > 



and < Simp < ln2. 



The anomalous dimension of the 



5 



auxiliary impurity fermionic field (see also App. is 



Vf 



3£ 



3£ 



3£" 



{i2 + '-E^)+<^i<^') for ,90^0 



for 50 = 



(14) 

r]f describes the decay of the two-point correlator of the 
/ fermions. In the Kondo model ^ this propagator is 
not directly measurable; however, it has been argued2& 
that mobile fermionic quasiparticles in a critical antifer- 
romagnet are described by a very similar theory, provided 
that their momentum is close to a minimum or van-Hove 
point of the band structure, such that the dispersion is 
quadratic and can be neglected. Then, rjf is the exponent 
of the electron Green's function as measured in photo- 
emission: 



for 7o 7^ if As = 0. The flow diagram. Fig. ^ shows 
that F-CR can be understood as a multicritical fixed 
point; consequently, universal one-parameter scaling in 
its vicinity is only realized for 70 — 0. For r — > 0, the 
F-CR fixed point merges with LM. 

NRG calculationsii have established that the fixed- 
point structure of the pseudogap Kondo problem changes 
for larger r, and two different critical fixed points occur 
for r > r* w 0.375. The particle-hole symmetric critical 
fixed point, present at small r, disappears for r > i. No- 
tably, the particle-hole symmetric and asymmetric crit- 
ical fixed points share many common propertiesiiiiSii^. 
In any case, the physics at large r is inaccessible within 
the present weak-coupling expansion; we will restrict our 
attention to the range of < r < ^ in the following. 
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(15) 



where cq is the (renormalized) energy of the quasiparti- 
cle. This situation, which can be regarded as the bosonic 
analogue of the familiar x-ray edge problem^^, applies to 
electrons in Kondo insulators at a magnetic critical point, 
and could also be relevant for holes near the anti-nodal 
points in high-temperature superconductors^^. 

As noted earlieri^iiS the B-FL fixed point is unstable 
w.r.t. breaking of the spin rotation symmetry. Further, 
for e ^ it merges with LM. 



5. Fermi-Bose critical fixed point (FB-CR) 

Finally, the FB-CR fixed point controls the transition 
between F-SC and B-FL, i.e., between the phase with 
fermionic screening and the bosonic fluctuating phase. 
For interacting bosons, ^ Q, we have 



*2 e 
^ = 2 



111 57r2\ er ^3 

1- Oie.rY , 

968 132 / 2 2 ^ ' ^ ' 



5 hn^ > ,0 

\ + OU.ry 

484 132 ' ^ ' ' 



(19) 



4- Fermionic critical fixed point (F-CR) 

The pure fermionic problem with a pseudogap, r > 0, 
features a critical fixed point with 



r=r + -+0(r3), 7* 







(16) 



which controls the transition between the decoupled and 
the Kondo-screened phases. As detailed below, we have 
< Tximp < 1/4 and 5'inip > ln2. The anomalous field 
dimension is given by 



3r2 



(17) 



As the fixed point is infrared unstable, we can discuss the 
scaling dimension of the leading relevant operator, which 
is just the inverse correlation length exponent of the LM 
- F-SC transition. Expanding the beta function around 
the fixed point value yields: 



(18) 



This exponent describes the vanishing of the character- 
istic energy scale T* when the system is tuned through 
the transition, i.e., T* cx |jo — iocV ■ 

The F-CR fixed point is unstable w.r.t a finite cou- 
pling to gapless bosons, therefore it cannot be reached 



whereas for 5 = we find 



J 



9 9 

6 e er r , o 

2 + ¥-y-y+^(^'^) ' 

r + |+C'(e,r)^ 



(20) 



The fixed point is again characterized by non-trivial uni- 
versal values of Tximp and S'imp depending on r and e 
only, in particular for r = we have < S'imp < ln2. 
The anomalous field dimension is 



^7/ 



f -¥(2i2+^)+^^(^'0' for.go^O 



0(e,r) 



for 50 = 



(21) 



The inverse correlation length exponent of the B-FL 
F-SC transition is 



1 



M363 + 132 j ~ -9 rr tor 50 7^ 



4£r 
9 



19r^ 
27 



for 50 = 
(22) 

up to terms of order Oit.rf' . Note that this exponent 
describes the behavior if the impurity model is tuned off 
criticality by varying jo or 70, while maintaining criti- 
cality of the bulk magnet (!), Ag = 0. If, in contrast, 
the bosonic bath is tuned away from its critical point, 
then the physics is dominated by Ag |SJ|, as discussed 
in Ref.Ei(see also Sec. HiTdI below). The FB-CR fixed 
point only exists for e > 0. 
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The results of Refs. ITlllT^ suggest that the fixed point 
character of FB-CR changes for larger r similar to that of 
the F-CR fixed point. In contrast, the properties of the 
bosonic part are likely smooth as function of e for < e < 
2 (1 < d < 3), as the (1 + e) expansion of Ref. yields 
universal properties similar to the (3 — e) expansion^; 
this is further supported by the large- theory presented 
in Ref. and App. |D| 



C. Phase diagram 

The schematic RG flow from this small (e,r) analy- 
sis is shown in Fig. The main conclusion is that the 
fermionic and bosonic Kondo couplings compete, i.e., the 
coupling of an impurity to collective spin fluctuations 
suppressed fermionic Kondo screening. This can be an 
important aspect for impurity moments in strongly cor- 
related materials-—. 

We expect the general structure of the phase diagram 
to be valid for all values < e < 2 and r > 0. This is sug- 
gested by numerics for the fermionic partifi^ii and various 
large- methodsi^iSS; for r > ^ particle-hole asymmetry 
is needed to reach the F-SC fixed point as noted above. 
For the case of a metallic fermionic bath, r = 0, the F- 
CR fixed point is absent, and the phase t rans ition line 
originates in the 7 = J = point, see Refs. I4l5l 



D. Tuning bulk vs. boundary criticality 

So far the analysis was restricted to s = Sc, i-e., the 
bulk quantum critical point of the host magnet; here we 
have seen that tuning the couplings jo and 70 can lead 
to boundary quantum phase transitions. 

It is instructive to also discuss the physics away from 
the bulk critical point, and in particular the behavior 
of the impurity properties upon tuning the bulk mag- 
net through its ordering transition. The behavior will of 
course depend on whether the values of the impurity cou- 
plings, jo and 7o, place the impurity problem at s = Sc 
into the F-SC or B-FL phases, or on the phase transition 
line. Fig. n 



1. r = 

For the metallic case, r = 0, a finite host spin gap, 
corresponding to s > Sc, always leads to an impurity 
screened by fermionic quasiparticles. 

Then, for large jo there is no qualitative change in 
the impurity properties when tuning (s — Sc) through 
zero, because for s < Sc the physics is described by a 
screened Kondo impurity in a weak magnetic field^. In 
contrast, for small jo the low-energy behavior of the im- 
purity changes from F-SC for s > to B-FL at s = s^j 
and for s < Sc the moment is polarized by the exchange 
field, with an effective magnetization cx (sc — s)'^*''^/^ 



(Ref. Il9j) . Only for jo , 70 placing the impurity problem 
on the transition line in Fig. ^ the impurity is driven 
critical by driving the bulk to the critical point. 

2. r > 

In the pseudogap situation, r > 0, the presence of a 
finite host spin gap, Ag > 0, still allows a transition 
between LM and F-SC. This can be understood in a two- 
step RG procedure: For UV cutoff energies larger than 
As the above RG can be apphed, generating a flow of j. 
Once the UV cutoff reaches Ag , the RG flow of 7 is ter- 
minated. The remaining RG procedure then effectively 
treats a model with fermionic bath only, but the flow of j 
starts at j{As) generated in the first part of the RG. Thus 
for r > there exists a transition between LM and F-SC 
as function of jo, controlled by the F-CR fixed point. 

We can now discuss the behavior when tuning the bulk 
magnet through its critical point. For large jo the im- 
purity is screened, and its properties remain essentially 
unchanged upon tuning (s — Sc) through zero as above. 
For very small jo the low-energy behavior of the impu- 
rity changes now from LM (!) for s > Sc to B-FL at 
s ~ Sc, and at s < the moment is polarized. Interest- 
ingly, there is now a significant intermediate range of jo 
values where the impurity is screened for large Ag, then 
becomes uncreened for smaller Ag (Ref. I20I) - this tran- 
sition is the one controlled by F-CR - and at s < Sc the 
behavior again is determined by B-FL^^. 

Summarizing this discussion, for the single-impurity 
model it is clear that, for generic values of the impurity 
couplings jo and 70, tuning the bulk magnet critical does 
not drive the impurity problem critical; coincidence of 
the two critical points requires additional fine tuning of 
parameters. 

IV. MAGNETIC SUSCEPTIBILITY 

The linear response functions to an applied field have 
been discussed for the purely bosonic model in Ref. ITgt 
here we supplement this by calculating the fermionic con- 
tributions and collecting the results for the different fixed 
points. Also, we briefly discuss hyperscaling properties 
with respect to a local field near the critical fixed points. 

The external field is coupled to the bulk magnetic de- 
grees of freedom using the following replacement in Ti^, 

TTq -> {TTa - iea/3.yHu/3ix)(l)^)'^ , (23) 

while we simply add field terms to Tic and Tiimp, 
-Hua{y){cl(7'^^Cf,){y) and 

^imp, a^a ■ (^^) 

The bulk field varies slowly as function of the space 
coordinate, and i?i,np is the magnetic field at the location 
of the impurity. 
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With these definitions, a spatially uniform field applied 
to the whole system corresponds to — i?imp = H. 
Response functions can be defined from second deriva- 
tives of the thermodynamic potential, — —TlnZ, in 
the standard way^«: x^^^ measures the bulk response to 
a field applied to the bulk, Ximp.imp is the impurity re- 



sponse to a field applied to the impurity, and Xu,imp is ity at criticality obeys 



B. Hyperscaling 

The structure of the RG implies that the impurity cor- 
relations obey certain hyperscaling properties at the non- 
trivial fixed points, including u/T scaling in dynamical 
quantities^. For instance, the dynamic local susceptibil- 



the cross-response of the bulk to an impurity field. Then 
the impurity contribution to the total susceptibility is 
defined as 



Ximp(T) = X- 



1, (25) 



imp, imp 



where Xu^u*^ is the susceptibility of the bulk system in 
absence of impurities. The local impurity susceptibility 
is given by 



Xloc(T) = Ximp,: 



(26) 



which is equivalent to the zero-frequency impurity spin 
autocorrelation function. 



A. Local susceptibility 

The local susceptibility Xioc(^^) at an intermediate- 
coupling fixed point is expected to follow a power law, 

1 



^un xioc(r) oc 



■Vx 



(27) 



This relation defines an anomalous exponent, 77^, which 
controls the anomalous decay of the two-point correla- 
tions of the impurity spin. (In Refs. [19II21I this exponent 
was labelled 77'.) 

Following the standard scheme^^, rj^ is calculated by 
introducing a xioc renormalization factor, Z^, evaluat- 
ing the lowest-order diagrams for xioc, and demanding 
the resulting renormalized expression to be free of poles. 
From the two-loop expression for Z^, given in Add. 1X1 



we find r]^ = 2(7^ - 
fixed points yields 

F-CR: = 
FB-CR: ??x = 

B-FL: = 



7^) Evaluating this at the RG 



(28) 



where the last two lines are for interacting bosons. We 
see no reason that the i]^ values for the FB-CR and B- 
FL fixed points should be equal to all orders. In contrast, 
for non-interacting bosons, g = 0, it has been shown in 
Rcfs. 5 19 that all higher-order corrections to 77^ vanish 
identically at a fixed point with finite 7*, thus 



FB-CR 
B-FL 



Vx 



(29) 



to all orders in perturbation theory. To the order we have 
worked here we find r/j = 3r/^/8, but this does not hold 
to higher orders. 



XL(-,T).-^<I>,(^) 



(30) 



where $1 is a universal crossover function (for the par- 
ticular fixed point), and Bi is a non- universal prefactor. 
The asymptotic behavior of $i(a;) is 



$i(x) 




(31) 



A similar scaling form applies for the T matrix, see below. 

In addition, for the critical fixed points F-CR and FB- 
CR we can also consider scaling properties away from, 
but close to, criticality. For details in the context of the 
pseudogap Kondo model we refer the reader to Ref. 0. 
In particular, the static local susceptibility follows 



xioc(r) 



B2 



$2 



jo - ioc 



(32) 



Using hyperscaling, the results for the correlation length 
exponent Sec. IIIII and for the anomalous exponent rj^ 
from above are sufficient to determine all critical expo- 
nents associated with a local magnetic fieldiS. In partic- 
ular, the T — > local susceptibility away from criticality 
obeys: 

Xloc(io > joc) OC (jo - joc)"^ , 7 = (1 - %) , 

Txioc(jo <joc) oc (joc - joV , j' = I'Vx- (33) 

With these definitions, the exponents 7 and 7'/2 coincide 
with the exponents 7 and /3 defined in Ref. 0. The 
perturbative expressions for the exponents in Eq. H33|) 
are 



F-CR: 7 = - + l+0{r) 
r 2 

7' = r + 0{r^), 



(34) 



FB-CR: 7 = 

7"' = 



2r + e 

2e 
2r + e 



+ 0(1), 
+ 0{e,r). 



At this order the results for g = and g ^ are identi- 
cal; the next-to-leading order contributions at the FB-CR 
fixed point can also be found from H22() . The results for 
the F-CR fixed point can be compared with numerical 
resultsi^, with excellent agreement, e.g., j{r = 0.1) = 
10.63 ± 0.02 from NRG, see also Fig. El 

As proposed earlieniiii^, limT^o 2^Xioc can serve as an 
order parameter, as it vanishes continuously as jo — > ioc 
from below, and is zero for jg > Joci i-e., in the F-SC 
phase. 
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a) 



20 



15 



10 




0.1 0.2 0.3 0.4 
r 



FIG. 2: Comparison of numerical data for the critical point of 
the pseudogap Kondo model (F-CR), obtained with NRG, to 
the analytical results of the renormalized perturbation theory, 
a) Correlation length exponent v and susceptibility exponent 
7 The analytical results are in (1181 and 13411 : to two- 

loop order, both exponents are equal, and sh own by the solid 
line. The symbols are NRG data from Refs. llOllSl b) Curie 
prefactor Cimp of the impurity susceptibility (1351 . Solid: an- 
alytical result 1391 : squares: NRG data taken from Ref. I 111 



C. Impurity susceptibility 

The quantity Ximp measures the impurity contribution 
to the total uniform susceptibility. Importantly, it does 
not acquire an anomalous dimension at an intermediate- 
coupling fixed point^, because it is a response function 
associated to the conserved quantity S'tot- Thus we ex- 
pect a Curie law 



lim Ximp 



(T) 



C'ie 



T 



(35) 



where the prefactor is in general a non-trivial univer- 
sal constant different from the free-impurity value 5(5-1- 
l)/3. Apparently, Eq. can be interpreted as the 

Curie response of a fractional effective spini^. 

Since the present expansion is done around the free- 
impurity fixed point of the 5 = ^ impurity, the result 
will always have the form Tximp = i + ^(T^Ximp)- Tech- 
nically, the absence of an anomalous dimension implies 
that no poles in e or r will appear in the perturbative ex- 
pression for Ximp, i-fi-, all poles from the different contri- 
butions have to cancel. Importantly, such pole-free con- 
tributions can also arise from the denominator in Eq. ||SJ), 
see App. ^ 

At this point a reference to the two-channel Kondo 
problem is in order, where it is known that Ximp does 
display an anomalous power law. This does not contra- 
dict the above statement, however, in the two-channel 
model the prefactor Cimp of Eq. ((^ actually van- 
ishes due to a "compensation" effect-^^, and the lead- 
ing low-temperature behavior arises from corrections to 
scaling22i2^. 

Proceeding with our calculation, we first focus on the 
fermionic part. The analysis of xioc has shown that the 
first corrections to free-moment behavior arise at order 
jg. This is different for Ximp^ there is a single diagram 



proportional to jo which contributes to Xu,imp: 

■ 1 I, .ir cosh"' (fc/2r) 
Xu,imp = -Jo^ y dfc|fc| (36) 

Evaluating the integral in the limit of infinite UV cutoff 
and taking r gives: 



A(Tximp) 



(f) _ _J 



(37) 



where we have already expressed the result in terms of the 
renormalized coupling j, jo — j/i"'^, and used {T/ nY — 
l + 0[r). 

The bosonic contributions have been discussed in 
Ref. 0; here we only quote the result for complete- 
ness. Importantly, the impurity susceptibility is not a 
well-defined quantity for non-interacting bosons {g = 0) 
because the application of a uniform field would lead to 
negative bosonic mode energies - this yields an infrared 
singular magnetic response of the bulk system, and a sim- 
ilar infrared singularity of Tximp- For interacting bosons, 
the lowest-order term givesi^ 



A(Txi,„p) 



(b,„t) 



(38) 



where the Hartree-Fock renormalization of the mass (0 is 
required to remove the infrared divergence of the integral. 

With these perturbative results, we can finally evaluate 
the impurity susceptibility at the intermediate-coupling 
fixed points, to lowest non-trivial order. For the bosonic 
part, the order e contribution can also be identifiedi^. 



F-CR: Txi, 



FB-CR: Txi, 



B-FL: Tximp 



r)+0{r^).. 




(39) 



9e 

T 




40 4 i ^ ' 



Notably, our result for the F-CR fixed point is in excellent 
agreement with the NRG of Ref. 11, see Fig. [3 

For the bosonic case no numerical results for small e 
are available; calculations^ for e = 1 (d = 2) have found 
a value for Tximp which is very close to the free-moment 
value J. We believe that this coincidence is due to large 
corrections from higher orders of the e expansion; the fact 
that our result for the fermionic F-CR fixed point agrees 
well with numerics strongly supports that the renormal- 
ized perturbation expansion employed here is the correct 
method to evaluate the impurity susceptibility (which 
had been questioned in Ref. 



V. IMPURITY ENTROPY 

In general, zero-temperature critical points in quan- 
tum impurity models can show a finite residual entropy 
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FIG. 3: Feynman diagrams for the impurity partition func- 
tion of the physical sector of the Hilbert space, i.e., {Q)x. 
These diagrams directly enter the calculation of the impurity 
entropy. Full/dashed/wiggly lines denote f/c/<j> propagators, 
and the dots are the interaction vertices, a) Unperturbed im- 
purity part, b) and c) Corrections from the fermionic and 
bosonic baths. 



- this is in contrast to bulk quantum critical points where 
the entropy usually vanishes at T = 0, but features an 
unconventional power law, S(T) cx T^. 

For the models at hand, the impurity contribution to 
the low-temperature entropy is obtained by a perturba- 
tive evaluation of the thermodynamic potential and tak- 
ing the temperature derivative. In the present expansion 
the result has the form Simp{T = 0) = ln(25'-f 1) + AS'imp 
with S 

It is worth emphasizing that the aim of the e expansion 
for the entropy is not to capture singular corrections (as is 
the case for observables which develop anomalous power 
laws), but to determine the power series for the coefficient 
of an observable with regular behavior. As shown below, 
fully universal quantities can be obtained by correctly 
interpreting the result of the perturbative calculation. 
The situation is similar to the evaluation of the impurity 
susceptibility in Sec. IIVI where the coefficient Cimp in 
Eq. H35|l was obtained in an expansion in e and r. 

The diagrams necessary for the evaluation of the ther- 
modynamic potential are shown in Fig. |21 In the fol- 
lowing we list the corresponding results; some details of 
the evaluation are described in Add. lO- in particular it 
is important to note that certain disconnected diagrams 
appear in the expansion of the thermodynamic potential. 

The coupling to the fermions. Fig. Osl, yields 

The integrals can be performed in the limit of infinite UV 
cutoff, with the result 

^Sl2, = f fo tan^ (2- 2i-^'-)r(2 + 2r)C(l + 2r) 

where r(z) and Q{z) are the Gamma and Riemann Zeta 
functions. In the weak-coupling regime the above term 
goes to zero as T^"" for T ^ 0, thus the entropy at the 
LM fixed point is In 2. However, in the quantum critical 



region, the dimensionless combination ji^T^ approaches a 
universal value at the F-CR fixed point, and a universal 
correction to the impurity entropy arises. This is easily 
seen by replacing jo with the renormalized coupling j, 
defined by jo = Using {T / = 1 + r\n{T/^i) + 

0{r^) and expanding in r one finally obtains: 



AS" 



(f) _ 



37r2 In 2 



imp 



.9 

jr. 



(40) 



Notably, this correction to S'imp is suppressed by an 
additional factor of r. We have verified that the contribu- 
tions from the next-order diagram, Fig.|3|32, proportional 
to j^, also receive an additional factor of r; thus the term 
(|40|l is the only one to order r^ at the F-CR fixed point. 

The lowest-order bosonic contribution from Fig. ^1 
turns out to vanish identically for the 5 = case, i.e., for 
non-interacting bosons, because the contribution of this 
diagram to the thermodynamic potential is temperature- 
independent. For interacting bosons, we take into ac- 
count the T-dependent mass renormalization arising from 
g, and obtain using Q in a similar manner as above: 



as: 



(b2,i, 
imp 




(41) 



For non-interacting bosons, non-vanishing corrections to 
the entropy arise only at order 7'* from the diagrams 
shown in Figs.|3|:2,|3t3. With Eq. (EU one finds: 



AS', 



(b4) ^ 
imp 



-9t[^(3c2) + (3c3) - (3cl)^ 
3^2 , 
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-7 



(42) 



With these perturbative results, we can finally eval- 
uate the residual impurity entropy at the intermediate- 
coupling fixed points, to the order available. 



F-CR: S 



imp 



ln2( l + ^rM +0(r5). 



FB-CR: Si, 



B-FL : Si, 




The last two lines are for interacting bosonsSS. In princi- 
ple, the two-loop result for the impurity entropy can be 
obtained from the diagrams in Fig. fusing the methods 
of Ref. 23; we do not attempt this here. The entropies 
for the FB-CR and B-FL fixed points get modified for 
non-interacting bosons (g = 0): 



FB-CR: Sn 



B-FL : Si, 



= ln2-|^.^ + 0(.^r3), 



In 2- 



(44) 



In App. ^we present a large- iV calculatiorJ^ of the en- 
tropy at the B-FL fixed point, which also yields an answer 
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of the form S'imp = ln2 — 0{e^) for small e, but provides 
results^ for all < e < 2. 

All of the above expressions are consistent with 

oF-CR > , r, > oB-FL 
•^imp - ill ^ - -^imp ' 

cF-CR ^ oFB-CR ^ c-B-FL . fAr) 
imp — imp — imp ' V ^/ 

these are the conditions that the impurity entropy al- 
ways decreases in the course of the renormalization 
group flow2i. Apparently, in the e expansion we have 



- S 



B-FL 



to the order calculated here, be- 



cJ^B-CR 
imp '"'imp 

cause the coupling constants 7 of both fixed points are 
equal to one-loop order. However, it can be expected 



that in general 5*; 



FB-CR 
imp 



B-FL 



due to the facts that 
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FB-CR 



< and ASl2p > 0- 



For the F-CR fixed point, the NRG calculations of 
Ref.O indicated a very small deviation of the impurity 
entropy from In 2. We have repeated these calculations 
and have found indications for a correction which scales 
as at small r, however, the accuracy is not sufficient to 
reliably determine the prefactor. Note that the entropy of 
the fermionic particle-hole symmetric critical fixed point 
also approaches^''- the value of ln2 for r ^ i, rendering 
the corrections to In 2 tiny for all < r < i. Interest- 
ingly, a similar entropy calculation can be done for the 
particle-hole asymmetric critical fixed point of the pseu- 
dogap Kondo model, where the expansion is performed 
around the upper-critical "dimension" r — 1. The cor- 



responding result is Si 
good agreement with numerics^ 



imp 



In 3 - (8/9)(l-r)ln2, 



VI. T MATRIX 

An important quantity in a Kondo model is the con- 
duction electron T matrix, describing the scattering of 
the c electrons off the impurity. It is useful to de- 
fine a propagator, Gt, of the composite operator T^- = 
flfa'Ca', such that the T matrix is given by T{uj) = 
joGri^^)- As with the local susceptibility, we expect a 
power law behavior of the T matrix spectral density at 
the intermediate-coupling fixed points: 

T(uj) cx , (46) 

note that at tree level Gt oc uj^ in the present pseudogap 
problem. 

We now evaluate the lowest-order diagrams of the T 
propagator in an expansion in 7 and j, introduce a Gt 
renormalization factor, Zt, and determine Zt by min- 
imal subtraction of poles. The result, given in lAlSfl . 
yields the anomalous exponent as rjT — — 2j4-j^ + 2(7^ — 
1' 



Evaluating this at the RG fixed points we obtain 



F-CR: r]T 
FB-CR: rjT 

B-FL : r]T 



-2r, 
-2r. 



(47) 



{i2 + T^)+0{e') (go^O) 



(go = 0) 




b) 



c) 



<3> 







FIG. 4: Feynman diagrams for the Green's function Gt en- 
tering the conduction electron T matrix. Notation is as in 
Fig. 121 the double line is the full / propagator, the shaded 
area the full j vertex, and the open dots are the external 
sources, a) Bare T matrix, b) and c) Corrections from the 
fermionic and bosonic baths, d) Full T matrix. 



Remarkably, the two-loop corrections for F-CR and FB- 
CR have cancelled. As announced in Ref. 20, this can- 
cellation holds to all orders. The argument parallels the 
one for anomalous exponent of the local susceptibility of 
non-interacting hosons, given in Refs. 5 18. For fermions 
without a relevant self-interaction, all diagrams for Gt 
can be obtained from the first diagram with full j ver- 
tices and full f electron propagators. This gives the rela- 
tion Z^^ = Zj /Zj; note that the corresponding relation 



for the local susceptibility of non-interacting bosons is 
Z-^ = Z?^/Zj. Using the definition of Zj (EHJ and tak- 
ing the derivative w.r.t. (In fi) at fixed bare Jq, we obtain 
= — r + — rjT 12. Thus, at any RG fixed point 

with finite j* we find the exact result 



^2r 



T(u;) 



cx UJ 



(48) 



which applies to the F-CR and FB-CR fixed points. For 
F-CR it is in perfect agreement with NRG result sA-'-- . 

At finite temperatures, hyperscaling again implies a 
scaling form of the T matrix, which can be written as: 



Bt 

T{u;,T) = ^^T 



(f) 



(49) 



vahd at both the F-CR and FB-CR fixed points. 

Away from criticality, scaling is also obeyed as function 
of the distance from the critical point, similar to Eq. H32|l . 
We note that in the LM regime, bare perturbation theory 
is applicable and gives T{lj) cx lj"^ . Then, for ^ joc in 
the pseudogap Kondo problem, one finds a crossover from 
T{uj) cx UJ-'' for w > T* to T(w) cx uo'' for w < T*, in 
agreement with NRG resultsiS. 

For transport experiments in a pseudogap system with 
dilute magnetic impurities, we can estimate the impurity- 
induced resistivity Pimp by evaluating the conduction 
electron bubble with the lowest-order self energy. The 
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above implies a divergence pimp cx: 7^-'"(i-'') at criticality 
(with logarithmic corrections at both r = and r = 1), 
whereas pimp vanishes according to pimp oc away 
from criticality. Here, the (1 — r) factor in the exponent 
arises from the momentum integration (consistent with 
a resistivity being independent of the scattering rate at 
r = 1, Ref. Hi). Thus, for < r < 1 and sufficiently 
close to the critical point, one expects a characteristic 
peak in the temperature dependence of pimp, located at 
the crossover temperature T* . 



VII. DYNAMICAL MEAN-FIELD THEORY 

Magnetic ordering transitions in certain heavy-fermion 
compounds, like CeCu6-a;Auxand YbRh2Si2, have been 
proposed to be accompanied by the breakdown of Kondo 
screening at criticality^''*'^. One particular scenario*' em- 
ploys an extended DMFT description for this transition, 
where the original Kondo lattice model is mapped onto 
a self-consistent version of the Fermi-Bose Kondo model, 
and the critical point of the bulk system is mapped onto 
the BF-CR critical point of the impurity model. Recent 
numerical calculations have found the existence of such 
a transition in the extended DMFT for an anisotropic 
Kondo lattice model^ (however, competing claims have 
also been published^). 

In the following, we want to focus on a property of 
such a locally-critical point which has not received much 
attention so far, namely the entropy and specific heat 
near the critical point. In principle, our statements will 
apply to all quantum critical points in DMFT where the 
effective impurity model is itself critical at the bulk tran- 
sition. (Note that this does not seem to be the case for 
the much-studied Mott metal-insulator transition in the 
one-band Hubbard model*.) 

To be specific, let us consider the thermodynamic po- 
tential of an Anderson lattice model, treated within ex- 
tended DMFT (equations for a t — J model are similar). 
It has been shown in Ref. that the lattice potential 
(per site) can be written as the sum of impurity and 
bath contributions, 

^ = r^hnp + + r?^ , (50) 

= T^^^ln[Gc(fc,ic^„)/Gioc(ic<^n)]e'""°^, 

^ k 

where Gc{k,iLOn) and D^(k,ivn) are the fermionic and 
bosonic bath Green's functions, respectively, and Ns is 
the number of lattice sites. Gioc and xioc are the (local) 
impurity electron propagator and susceptibility - in our 
notation, Gioc is the propagator of the / fermions. The 
lattice entropy (per site) is simply — cJ^fJ. The temper- 
ature derivatives of the two terms f2c and vanish in 
the zero-temperature limit: note that they just represent 



contributions from systems of free fermions or bosons, 
where the respective densities of states are (at most) log- 
arithmically singular within the locally-critical point sce- 
nario in two space dimensions^. 

Thus, the entropy of the impurity model, —dT^irap^ is 
the only contribution to the total entropy of the lattice 
model which survives as T ^ 0. The results of Sec. [V] 
suggest that the critical point of the Fermi-Bose Kondo 
model, BF-CR, generically has a finite residual entropjiSi, 
with the perturbative result given in Eqs. 143144(1 . Ac- 
cording to (|50|1 . this implies an extensive entropy of 
the Kondo lattice model described by extended DMFT, 
which should be present at the T = critical point as well 
as in the low-temperature limit of the quantum critical 
region. As the ground state entropies of the stable phases 
(heavy Fermi liquid and antiferromagnet for the case of 
the Kondo lattice) are non-extensive, a huge anomaly in 
the specific heat, G = T{dS/dT), should be observed 
upon crossing the phase diagram at low temperatures 
above the quantum critical point. To our knowledge, this 
has not been observed experimentally. Furthermore, it 
can be considered unlikely that a quantum critical point 
of a lattice model displays an extensive residual entropy, 
as such a point would be unstable to any perturbation. 
(Note that a generic lattice quantum critical point has 
only one unstable direction in the RG flow, and vanish- 
ing residual entropy.) 

Finally, as detailed in Ref. 20, the DMFT self- 
consistency equation for the fermionic bath at local crit- 
icality can only be fulfilled with a fermionic spectrum 
with r = 0. Therefore, the local critical point displays 
at most logarithmic corrections to a metallic local den- 
sity of states, and one can thus expect a logarithmic T- 
dependence of the resistivity. A discussion of the resis- 
tivity away from criticality requires to consider the fully 
self-consistent problem, and is beyond the scope of this 
paper. 



VIII. DISCUSSION AND CONCLUSIONS 

We have calculated critical properties of a Fermi-Bose 
Kondo model, describing an impurity spin coupled to 
both fermionic quasiparticles with a pseudogap DOS 
and collective mode bosons, using renormalization group 
techniques. In particular, we have determined the im- 
purity susceptibilities and residual entropies at different 
intermediate-coupling fixed points. 

Our results for the phase transition in the pseudogap 
Kondo model can be compared with numerical results 
obtained by high- accuracy NRG calculations. The excel- 
lent agreement strongly suggests that the renormalized 
perturbation expansion approach for critical properties is 
reliable in the context of impurity quantum phase tran- 
sitions (provided that the fixed point is in the applicabil- 
ity range of the renormalization group procedure), and 
gives numerically accurate results for a finite range of the 
(small) parameters r, e. 
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Remarkably, the physics found in the present models 
is somewhat different from the one of the two-channel 
Kondo effect. In the two-channel problem, a num- 
ber of thermodynamic properties vanish in the strict 
scaling limit, and corrections to scaling yield the lead- 
ing low-temperature observables. In contrast, for the 
intermediate-coupling fixed points considered here, the 
predictions of naive scaling apply. The reason for this 
difference apparently is the finite tree-level scaling di- 
mension of the impurity couplings in the present models, 
whereas the coupling in the metallic two-channel problem 
is marginal. 

As can be seen from the perturbative results presented 
in this paper, the prefactors of most of the higher-order 
corrections in the e expansion of physical quantities are 
not small, and thus no quantitatively accurate numbers 
can be obtained for e = 1. It may be possible to use 
certain resummation techniques of the perturbative ex- 
pansion - this is beyond the scope of this work. Never- 
theless, we expect that the qualitative physics is captured 
by the present study, which is supported by the large- 
calculations in Ref. .19. and App.^for the purely bosonic 
Kondo problem. 

On the single-impurity application side, the mod- 
els discussed here can describe magnetic impurities in 
strongly correlated materials, where the interplay of 
quasiparticles and collective low-energy excitations is a 
central ingredient. As has been argued elsewhera^2i2S44, 
the pseudogap Kondo model appears relevant for impu- 
rity moments in high-temperature superconductors. De- 
tailed calculations of the local density of states in the 
vicinity of the impurity^^iM have been compared with 
scanning tunneling microscopy results^ in Zn-doped 
Bi2Sr2CaCu208+i. In addition, the interplay with bulk 
spin fluctuations possibly explains the large doping de- 
pendence of the Kondo temperature across the phase 
diagram of YBa2Cu3 07_,fi^4i^, as has been explicitely 
demonstrated in calculations using the pseudogap Fermi- 
Bose Kondo model in Ref. 20. In this case, the inclusion 
of the bosonic self- interaction is important, as fermionic 
quasiparticles and spin fluctuations are decoupled in the 
low-energy physics of a d-wave superconductor, and Lan- 
dau damping is absentia. 

Finally we have discussed the issue of the ground- 
state entropy in scenarios of local quantum criticality, 
described by dynamical mean-field theory. We found 
that the finite residual entropy of the critical impurity 
model, found in Sects'", leads to an extensive ground- 
state entropy of the lattice model, which appears unlikely 
to be realized in a finite-dimensional system. One can 
conclude that the dynamics at a heavy-fermion critical 
point, which is associated with the breakdown of Kondo 
screening, cannot be entirely local, as described by dy- 
namical mean-field approaches. In other words, at lowest 
temperatures the entropy has to be quenched by non- 
local correlations. Versions of DMFT can still provide 
a reasonable description of the physics at intermediate 
energies and temperatures. 
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APPENDIX A: FIELD-THEORETIC RG TO TWO 
LOOPS 

The renormalization group to two-loop order is con- 
veniently done in the field-theoretic scheme, employing 
dimensional regularization and minimal subtraction of 
poles. Using the conventions of Ref. |2^ we introduce 
renormalized fields and dimensionless couplings for the 
bulk bosons 



50 



= V Z (j)Ra 



-9 



and for the impurity 

70 = 
jo -- 



i'/^Z^ 



ZSd+i 



^z. 



Zf 



-3 ; 



(Al) 
(A2) 

(A3) 
(A4) 

(A5) 



here /i is a renormalization energy scale, and 
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Sd 
Sd 



r(d/2)(47r)'*/2 

r(d/2- 1) 



^^d/2 



(A6) 



No renormalizations are needed for the fermions as their 
self-interaction is assumed to be irrelevant in the RG 
sense. 

The needed diagrams have been evaluated in Refs. 
and|^ and will not be repeated here. Minimal subtrac- 
tion of poles yields the renormalization factors for the 
bulk bosons'^^ 



Z 



= 1 - 



144e ' 

11.9 /121 
6e V 36e2 



37 
36e 



(AT) 



and for the impurity 
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The RG beta functions can now be evaluated by taking 
the n derivatives of HA2|I . (|A4() . and (|A5() at fixed values 
of bare coupling constants. For the bulk coupling we find 
the known result: 



Kg) = M 



d5 
d/i 



= -eg 



llg2 23.g3 
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which has the stable fixed point 

414e2 



5 = TT 



6e 
11 



1331 



(A9) 
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The impurity couplings are governed by the following 
beta functions: 
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-J7 -J7 



36 



■ 9 

J 7 



(All) 



Using these beta functions one obtains the fixed point 
values of the couplings quoted in Sec. IIII Bl 

The anomalous field dimension ?y/ is obtained from 



rjf ^fi 



din Zf 



(A12) 



where the derivative is evaluated at the RG fixed point. 

To calculate the anomalous exponent of the local sus- 
ceptibility, a renormalization factor, Z^, for the two- 
point correlations of the impurity spin has to be intro- 
duced. The required diagrams have been evaluated in 
Refs.iandlii the result for Z^^ within the minimal sub- 
traction scheme is 



Z^^l- 



27^ 



7_ 
e 



3_ 

2r 



Employing 
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dfj, 



(A13) 



(A14) 



one obtains the rj^ values quoted in Sec. IIVI 

Similarly, the renormalization factor for the T matrix, 
discussed in Sec. I VII is obtained as 



Zj 
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2r 
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e 
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APPENDIX B: ROLE OF (Q) 

Here we briefly discuss the role of the denominator 
appearing in the equation for observables, ((HJ, within 
the pseudo-fermion technique employed here. (Note that 
a similar denominator also appears in the formalism used 
in Ref. 19.) 

For a renormalization group procedure which directly 
calculates renormalizations of vertices and propagators 
(as the momentum shell method), the denominator does 
formally not appear. In contrast, when calculating true 
observables, as e.g. done in Ref. 0, (Q) needs to be 
taken into account. Importantly, approximating the de- 
nominator with its leading contribution, 2exp(— /3A), is 
in general not sufficient, as the denominator contains - 
as the numerator of ||HJ) - contributions from all orders in 
an expansion in the non-linear couplings. 

However, in the calculation of critical exponents using, 
e.g., the minimal subtraction scheme, one only needs to 
keep track of contributions which develop poles in e after 
dimensional regularization, i.e., are logarithmic at the 
marginal dimension. Typically, the denominator does 
not develop such poles, which can be seen by compar- 
ing power counting for response functions, i.e., second 
derivatives of the thermodynamic potential, and for the 
potential itself. Therefore, for the calculation of anoma- 
lous exponents it is permissible to ignore the perturbation 
expansion of (Q) , as done in Ref. |a. 

In contrast, when calculating a quantity which does 
not develop an anomalous power law, like Ximp in the 
present impurity problems, all poles cancel, and one has 
to collect contributions non-singular in e. Those also arise 
from a perturbative expansion of the denominator, and 
typically cancel corresponding terms from the numerator. 
Explicitly, we have to lowest order: 



(Q) = 2cxp(-A/3)(1-— 7^ 



d'^k 1 



(27r)" 



(Bl) 



to this order only the bosonic part enters, and eu is the 
boson dispersion appearing in D^{k,ujn) = [uj"^ + e^)^^. 

Finally, we note that in certain impurity problems the 
real part of the propagator renormalization requires the 
introduction of a suitable counter-term - physically this 
reflects the shift of the phase transition point, arising 
from the real part of a self-energy. In the present case no 
counter-terms are necessary. 



APPENDIX C: PERTURBATIVE EXPANSION 
FOR THE IMPURITY ENTROPY 

In this appendix we describe the perturbation theory 
for the impurity part of the thermodynamic potential. 
The impurity spin is represented by / pseudo-fermions as 
above, together with the chemical potential in the limit 
A ^ oo. We shall illustrate that the correct treatment 
of this limit has non-trivial consequences for the higher- 
order diagrammatic expansion of the entropy. 
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The starting point is the partition function in the phys- 
ical sector of the Hilbert space (Q = flfu = 1), which 
can be written as 



^imp 
^imp,0 



hm 



Ml 

{Q) 



(CI) 



A,0 



Here, (. . .)a is an expectation value with the full Hamil- 
tonian in the presence of A as above, and (. . .)x,o is the 
expectation value in the absence of the coupling between 
impurity and baths. The above expression is easily un- 
derstood: for A oo, {Q)\ is the partition in the physi- 
cal sector times exp(— A/3), i.e., the "factor" Q suppresses 
the = part of the Hilbert space. Further, in this limit 
(Q)a.o = 2 exp(— A/3), and the unperturbed impurity part 
of the partition function is just .^imp.o = 2 for a spin-^ im- 
purity. Importantly, the limit A ^ oo in Zi^p suppresses 
the contributions from disconnected diagrams, consist- 
ing of more than one loop of / pseudo-fermions, because 
those are smaller by further powers of exp(— A/3). 
The thermodynamic potential is given by 



^imp,o — —Tin ■ 



^iinp,0 



(C2) 



The partition function IjCip can now be expanded in the 
impurity couplings using standard diagrammatics, and 
then the log in (|C2|I is expanded to the required order. 
There is now no cancellation of disconnected diagrams in 
^^imp, because those diagrams do not appear in Zi^ip ~ 
this reflects the fact that we are not expanding around a 
system of free particles here. 

For the bosonic impurity problem, we have 



^imp 
^imp.O 



1 + 



70^ 



^(4) 



(C3) 



where is the sum of all n-th order different connected 
diagrams; 

i:>(2) and are shown in Fig.[3|:;. Thus, 



^IlM ^ 20 n(2) , To n(4) „ lQ_r){2f , 

-T 2 4 8 ' 

(C4) 

i.e., the expansion for the thermodynamic potential con- 
tains disconnected diagrams! 

The entropy is obtained from fiimp by Simp — 
— 9rJ7ii„p. Power counting shows that ^imp = In 2 as 
T ^ in the cases where the fixed point value of the 
impurity coupling is zero (as is the case, e.g., in the 
Bose-Kondo model for d > 3 space dimensions). Then, 
only subleading corrections arise from the coupling to 
the bath. In contrast, in the case of a finite renormalized 
coupling j* or 7*, bare perturbation theory is infrared 
divergent, as for standard problems below their upper- 
critical dimension. Proceeding in the scheme of renor- 
malized perturbation theory, the perturbative result (to 
lowest non-trivial order) can be expressed in terms of the 
renormalized couplings. At the RG fixed point, the cou- 
plings can then be replaced by their fixed-point values, 
and the result for the entropy is finite as T ^ 0. 



The impurity part of the thermodynamic potential di- 
verges with the cutoff, i.e., we have riimp = -Eimp — T^imp, 
where -Bimp is the non-universal (cutoff-dependent) im- 
purity contribution to the ground-state energy. However, 
the impurity entropy 5i,np is fully universal, and the UV 
cutoff can be sent to infinity after taking the temperature 
derivative of r2i,nn. 



APPENDIX D: LARGE-TV THEORY OF THE 
BOSE KONDO PROBLEM - ENTROPY 

As the perturbative method employed in the body of 
the paper is restricted to small e = 3 — d, alternative 
approaches are desirable which can provide results for all 
e. In Ref.0 a dynamic large- approach to the bosonic 
Kondo problem has been developed. Here we shall use 
this approach to evaluate the impurity entropy at the 
B-FL fixed point in arbitrary space dimension d, which 
complements the results of Sec. |V] 

We begin by summarizing the large-A^ theory of 
Ref. [19I The impurity spin symmetry is generalized 
from SU(2) to SU(A^). The spin operator is represented 
by auxiliary fermions fu {i^ — Ij-.-iA^), and a chemi- 
cal potential Aq is introduced to enforce the constraint 
flfu = NqQ. We will restrict the considerations to 
go — 1/2; this choice has the advantage of preserving 
particle-hole symmetry, such that the chemical potential 
Ao is actually zero for all temperatures. The Hamilto- 
nian is given by the natural large-A" generalization of 
'^imp + "^0 of Sec. m where the (f) field now has A^^ — 1 
components; for details see Ref. ^3 Taking the limit 
N ^ 00 results in a dynamic saddle point^, and the 
impurity physics is captured by a self-consistent integral 
equation. 



^f[T)^-^lD^{T)Gs{T), 



(Dl) 



where is the local spin fiuctuation propagator lO, 
and the self-energy S/ is defined by: 



Gf^iiujn) = it^n - Ao - S/(iw„) 



(D2) 



Here Gf{T) — — (T/(r)/^^(0)) is the /itZ/ auxiliary fermion 
Green's function in standard notation. Clearly, IjDll) 
represents a summation of all self-energy diagrams with 
non-crossing boson lines, somewhat similar to the non- 
crossing approximation known from fermionic Kondo 
physic o^'^- i. In writing down the above equations we have 
assumed full SU(A^) symmetry and dropped the corre- 
sponding SU(A'^) indices. Finally, Aq is determined by 
the constraint equation: 



G/(t = 0-) = i^G/(*c.„)e'^"""^ - Nqo, 



which is designed to satisfy the contraint f^fi, = Nqq on 
average. 
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FIG. 5: Large-A'^ result IDlOl for the residual impurity en- 
tropy (per flavor) at the B-FL fixed point, valid for a bath of 
non-interacting bosons in (3 — e) space dimensions. 



In the following we concentrate on the bulk critical 
point, s — Sc, where the bosonic spectrum is gapless, 
p^{uj) cx sgn(w)|cj|^~'^. It has been established in Ref. [l9l 
that the T = result for the fermionic Green's function 
obeys: 



Gf{r) 



1 



e/2-1 



(D3) 



in the long-time or low-energy limit; here Pf{uj) is the 
spectral density of G/. For the case of non- interacting 
bosons it is possible to obtain an analytical result for the 
finite-temperature scaling form of Gf, 



Gf{r) 



-A- 



ye/2 



7o \sin7r(TT) 



^/2 



(D4) 



which is dictated by conformal invariance^^ ; A is an am- 
plitude prefactor depending on e. 

The local susceptibility in the large- limit. 



Xlo 



r dTGf{T)Gf{-T), (D5) 
Jo 



can be evaluated directly at the critical point, s = Sc, 
with the result xioc oc r~^+'^. Thus we have the large- 
result 



Vx 



(D6) 



which coincides with the exact result H29|) for non- 
interacting bosons in the SU(2)-symmetric Bose Kondo 
problem. 

Now we proceed with the calculation of the residual 
impurity entropy, following the method by Parcollet et 
fflZ.^*^ for the fermionic NCA. At the saddle point, the 
impurity part of the free energy per flavor can be written 
as: 

/imp = go Ao + T ^ In Gf {iujn) - / dr S/ (r)G/ (-t) . 

(D7) 



For particle-hole symmetry, qo = 1/2, we have Aq = 0. 
Furthermore, the last term in (|D7|I behaves as const + 
0(T^+^/^), and thus does not contribute to the zero- 
temperature entropy. It remains to evaluate the con- 
tribution from Tr In G/. As shown in Ref. i46i. it can be 
written as: 



/imp - / dujufiu) 

TT 



arctan 



G}H 
G'Ucj) 



, (D8) 



with Gp Gy being the real and imaginary part of the 
finite-temperature retarded fermion Green's function, 
and n/ is the Fermi factor. 

In the following we restrict our attention to the case 
of non-interacting bosons which allows to make analyt- 
ical progress. The temperature derivative of the above 
expression has to be taken using the scaling form (|D4|I of 
Gf. After straightforward algebra in analogy to Sec. VI 
of Ref. 46 we find for the impurity entropy (per flavor) : 



^imp 



- — / ^ [u arctan(Mi) — arctan t] , (D9) 

n In 1 - U'^ 



where the substitutions l/t = tan(7re/4) and u = 
tanh(a;/2T) have been employed. The above result can 
be re-written as 



2 f^^\ uarctan(l/u) 
Si,„p = In 2 - - / dv ^ , V 



(DIO) 



This represents our large- iV result for the impurity en- 
tropy at the B-FL fixed point for non- interacting bosons. 
Numerical evaluation of the integral yields Simp as shown 
in Fig. 13 For small e, where l/t — 7re/4, the correction 
to the log 2 entropy is quadratic in e, which is similar to 
the SU(2) result in (gj). In particular, we see that simp 
is finite for all < e < 2, with the numerical value 



apid 



at e = 1. (Note that, e.g. 



0.4648477 



(Dll) 



only singular for < e < 1.) Sjmp 



the local susceptibility is 
vanishes as d — > 1, 
i.e., at the lower critical dimension of the bulk theory, 
where Gf is no longer singular - this is consistent with 
the considerations of Ref. |23. 

Unfortunately, the present large- iV theory cannot be 
easily generalized to the Fermi-Bose Kondo model with 
a single fermionic screening channel. Instead, the natu- 
ral dynamic large- iV formulation of the fermionic sector 
leads to a multichannel Kondo model^, which, however, 
does not display a strong-coupling fixed point with van- 
ishing residual entropy. 
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